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Abstract 

We study the spectrum of the QCD Dirac operator near zero virtuality. We argue 
that it can be described by a random matrix theory with the chiral structure of QCD. 
In the large limit, this model reduces to the low energy limit of the QCD partition 
function put forward by Leutwyler and Smilga. We conjecture that the microscopic limit 
of its spectral density is universal and reproduces that of QCD. Using random matrix 
methods we obtain its exact analytical expression. This result is compared to numerically 
calculated spectra for a liquid of instantons, and we find a very satisfactory agreement. 
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1. Introduction 

As has become clear from a wide body of botli numerical and theoretical evidence, the 
SU{Nf) X SU{Nf) chiral symmetry of the massless QCD Lagrangian is broken sponta- 
neously to SU{Nf), whereas the Ua{^) symmetry is broken explicitly by the anomaly (see 
The order parameter that characterizes the presence of chiral symmetry is the chiral 
condensate. According to the Banks-Casher formula , it is proportional to the spectral 
density of the Euclidean Dirac operator near zero virtuality. Therefore, the behaviour of 
the spectrum in this region is of fundamental importance to the phenomenon of chiral 
symmetry breaking. 

In this paper we study the spectrum of the QCD Dirac operator near zero virtuality in 
the thermodynamic limit while keeping the average spacing between the eigenvalues fixed, 
i.e. the spectrum is blown up by a factor V4 as V4 — > 00. This limit of the spectral density, 
which has a nontrivial limit as a consequence of a finite value of the chiral condensate, 
will be called the microscopic spectral density. It measures correlations on the order of 
one level spacing as opposed to the average spectral density that is obtained by averaging 
over many levels. We conjecture that the microscopic spectral density is universal. This 
is reminiscent to the universality observed in the theory of critical phenomena where large 
scale fluctuations wash out the details of the interaction. Because the spectrum near zero 
virtuality is dominated by the soft modes of the system, we expect that it is entirely 
determined by the symmetries of the microscopic theory. Indeed, Leutwyler and Smilga 
1^ were able to express negative moments of the microscopic spectral density in terms of 
the effective low-energy limit of QCD that only embodies its underlying chiral structure. 

A similar splitting of scales between spectral fluctuations on the order of one level 
spacing and the average spectral density is well-known in the statistical theory of spectra 
1^ and quantum chaos 0. In that case one finds that correlations on the microscopic 
scale are universal if the classical counterpart of the system is chaotic. The correlations 
only depend on the symmetries of the system and are given by the invariant random matrix 
ensembles. This suggests to describe the microscopic spectral density of the QCD Dirac 
operator in terms of a random matrix theory. The relevant ensemble will be called the 
chiral random matrix ensemble. We have shown that it is equivalent to the low-energy 
limit of the chiral effective Lagrangian (see section 3). Moreover, techniques developed in 
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random matrix theory enabled us P] to obtain an analytical expression for the microscopic 
spectral density (see section 4). 

Apart from the above arguments, there is additional evidence in favor of the univer- 
sality of the spectrum near zero virtuality. First, the Leutwyler-Smilga sum rules are 
not only found in QCD but in the Schwinger model as well and, numerically, for an 
instanton liquid model of the QCD- vacuum Second, the chiral random matrix model 
reproduces all Leutwyler-Smilga sum rules [0. Third, the microscopic spectral density 
of a Dirac operator in a liquid of instantons agrees within numerical accuracy with the 
random matrix result (see section 6). Fourth, the same microscopic spectral density is 
found in the Hofstadter model for universal conductance fluctuations [|10] . 



The Leutwyler-Smilga sum rules and the microscopic spectral density are based on 
the observation that a nonzero value of the chiral condensate implies that the spacing 
of the small eigenvalues is ~ I/V4, as opposed to ~ ^^/V^^^ for a noninteracting system. 
Such accumulation of small eigenvalues is obtained naturally for a liquid of instantons, 
where the exact fermionic zero modes can be chosen as basis states to construct a matrix 
representation of the Dirac operator |]lT|. Such model can be simulated numerically and 



results for its spectral density are presented in section 7. 

This paper is organized as follows. In section 2 we give an outline of the general 
framework of this paper and introduce the microscopic spectral density. A derivation of 
the simplest Leutwyler-Smilga sum rule is given in section 3. In section 4 we formulate 
the chiral random matrix model for the trivial topological sector of the QCD vacuum. 
Its equivalence to the low-energy effective partition function is shown in section 5. The 
microscopic spectral density is calculated in section 6, and, in section 7, numerical results 
for a liquid of instantons are presented. Concluding remarks are made in section 8. 



2. Formulation of the problem 

In this paper we study the spectrum of the massless Euclidean Dirac operator defined 
by the eigenvalue equation 
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(Z75 + 7A)0A = (2.1) 

where A is an SU (Nc) valued gauge field. The statistical distribution of the eigenvalues 
is induced through the gauge field configuration by the statistical weight of the Euclidean 
QCD partition function 

Z = Y.e-' < n n (An + m))m) >s^^a), (2.2) 

1^ /=1A„>0 

where the masses play the role of sources. The product is over the positive eigenvalues of 
the Dirac operator, and the average < ■ ■ ■ >s^{A) is over gauge field configurations with 
topological quantum number u weighted by the gauge field action S^{A). The topological 
part of the action, exp i6i>, has been displayed explicitly. The eigenvalues occur in pairs, 
A and —A, with corresponding eigenfunctions (px and 750a- The only exception is A = 0, 
when we have the logical possibility that 750a ~ 0a- Indeed, according to the Atiyah- 
Singer theorem, gauge field configurations with topological quantum number u have at 
least u zero modes with a definite chirality. 
The propagator is given by 

si.y)^!:-'-^. (2.3) 

where we have suppressed the color and Dirac indices. The chiral condensate < gg > is 
defined as the space-time average of < TrS'(x, x) > (the brackets < ■ • ■ > without indices 
denote averaging with respect to the partition function under consideration). According 
to the Banks-Casher 0] formula we have 

< gg >= -1 / d^x < TtS{x, x) >= i7f ^ ^^^^ ^ , (2.4) 

where the spectral density p(A) is defined as 

p{\)=Y.S{\-Xn). (2.5) 

An 

The condensate can be expressed as a derivative of the partition function 

i d 

<qq>= lim lim — logZ(m/). (2.6) 



The spectral density at zero virtuality has to be defined with care: the thermodynamic 
hmit has to be taken before the chiral hmit, i.e., 

1 TTl 

p(0) = hm hm - V -. (2.7) 

^ m^OVi^oo IT ^ + 

In this way it is possible to obtain a nonzero spectral density < p(0) > even in the case 
that zero eigenvalues are absent at all finite volumes. If this leads to a nonzero value 
of the chiral condensate, this phenomenon is referred to as the spontaneous breaking of 
chiral symmetry. A necessary requirement for this to happen is that 

< p(0) >~ \/4, (2.8) 

or, put differently, the spacing between the eigenvalues near zero virtuality is ~ I/V4. 

This number should be contrasted with the spacing of the small eigenvalues for a 
non-interacting system which is ~ (obtained from the solution of the free Dirac 

equation in a hypercubic box of volume V4). A natural explanation for the accumulation 
of eigenvalues near zero is provided by the properties of gauge field configurations of a 
liquid of instantons. Each isolated instanton, a localized field configuration by itself, gives 
rise to one localized zero mode. When the space-time volume is saturated with an equal 
number of instantons and anti-instantons the low-lying modes can be expressed as linear 
combinations of the exact zero modes and form the zero mode zone [0. For a rather 



dilute ensemble we therefore expect that the number of small eigenvalues is ~ V4. 

As was observed by Leutwyler and Smilga, eq. (2.4) implies the existence of a family 
of new sum rules. The simplest one involves the sum 

E < ^ >^ (2-9) 

which should converge to a finite limit for V4 — * cxo. 

The most important consequence of chiral symmetry breaking is the appearance of a 
massless pion in the chiral limit (m 0). In other words, the pion correlator at zero 
momentum (note that S{y,x) = — 755'"''(?/, x)75) 

n(p = 0) = J d%< TtS{x, y)S\x, y) > (2.10) 

diverges as ~ 1/m. This follows from the chiral Ward identity 

U{p = 0) = ^(< TiS{x, x)> - < TtS\x, x) >), (2.11) 



which allows us to express this correlator in the spectral density: 

n(p = 0) = ^-TT < p(0) >, (2.12) 
mv/4 

where we have used the translational invariance of < Tr5'^(x, x) >. From eq. (2.4) it 
then immediately follows that a finite value of the condensate leads to a massless pion in 
the chiral limit. As always, it is understood that the thermodynamic limit is to be taken 
before the limit m — 0. 

Because of the importance of the small eigenvalues in the spontaneous breaking of 
chiral symmetry, we propose to study the spectrum near zero virtuality in the thermo- 
dynamic by enlarging it proportional to V4. In this way the spacing of the eigenvalues 
remains 0{1) for V4 — 00. We call this limit the microscopic limit. Formally, it is defined 
as 

ps{x) = hm i <p{^) (2.13) 

in the sector of topological charge u. All Leutwyler-Smilga sum rules can be expressed in 
the microscopic spectral density and similarly defined spectral correlation functions. For 
the sum (2.9) we find 

J,.^. (2.14) 

By using random matrix theory we will obtain an analytical expression for ps{x) in the 
zero topological charge sector. 



3. Derivation of Leutwyler-Smilga sum rules 

According to Leutwyler and Smilga P|, the mass dependence of the low-energy limit 
of the QCD partition function is given by 



zim = 0, 



/ VUexp(^TT(mU-^exp(-ie/Nf) +mUexp(ie/Nf)). (3.1) 

JdetU=l 2 



The integration is over SU{Nf) with the Haar measure, m is the quark mass matrix, 
which is taken diagonal, and S is the the absolute value of chiral condensate. For one 



flavor, when the integration is absent, this ratio simphfies to 

J^'^'l^, = exp(m\/4Scos^). (3.2) 
Z[m = 0, 6) 

Expanding this ratio in powers of m and comparing the coefficients to those of the same 
expansion of the full QCD partition function provides us with an infinite family of sum 
rules. The simplest sum rule is derived from the 0{m?) term. In the zero topological 
charge sector (projected out by integrating over 0) the result is 



where the case Nf = 1 follows immediately from eq. (3.2). 



4. The chiral random matrix model 

In this section we construct a random matrix model for the spectrum of the Dirac 
operator near zero virtuality in the sector of zero total topological charge. The underlying 
assumption is that the microscopic spectral density is universal, and only follows from the 
symmetries of the system. For A'^^ flavors, the partition function that reflects the chiral 
structure of QCD is given by [0, 0, 0] 

Nf 

det ( '"^ ' 



Z = /pTP(T)ndet f , (4.1) 



where the integral is over the real and imaginary parts of the matrix elements of the 
arbitrary complex N/2 x N/2 matrix T, i.e., VT is the Haar measure. In agreement 
with the maximum entropy principle []T6[ the distribution function of the overlap matrix 
elements -P(T) is chosen Gaussian 

/\rv2 

P{T) = exp( —TiTT^). (4.2) 

The symplectic structure is a manifest consequence of chiral symmetry, and implies that 
the quark eigenvalues occur in pairs. The density of the modes that define the matrix 
representation of the Dirac operator, N/V^, is taken equal to 1, which allows us to identify 
the space-time volume and the total number of modes A^. 



The partition function (4.1) is inspired by the instanton hquid approximation to the 



QCD- vacuum [jT^, |T2[, where the fermion determinant is approximated by its value in 
the space of the fermionic zero modes. In that case, the overlap matrix elements Tij 
are a function of the collective coordinates of the instantons. Therefore, they are not 
statistically independent, and, as we will see in section 7, the average level density differs 
from a semi-circle, the result that can be obtained for the ensemble (4.1). 

Although in this section the random matrix model has been defined for zero topological 
charge only, it can be easily extended to the general case including an arbitrary value 
of the QCD angle. The essential observation is that if T is an x M rectangular 
matrix, the matrix in eq. (4.1) has exactly — M\ zero eigenvalues for ruf = 0. For 
simplicity we restrict ourselves in this paper to square overlap matrices. 



5. Reduction of the partition function 

In order to evaluate the partition function (4.1) the determinant is written as an 
integral over Grassmann variables 

( r/, 'I ) ^ /n^^wexpE ( ) ( - - ) ( ^ ) . 



/ 



(5.1) 



where the measure of the Grassmann integration is as usual 

Vijf = l[dij{di:t, (5.2) 



and the conjugation * is the conjugation of the second kind (i.e., ip** = [HI) Hll- The 
integral over T is gaussian and can be performed easily. In the partition function this 
results in the factor 

2 

which represents a 4-fermion interaction. 

The quartic term can be written as a sum of two squares 



(5.4) 
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Each of the two squares can be hnearized with the help of a Hubbard-Stratonovich trans- 



formation [T^. This allows us to perform the Grassmann integrations at the expense of 
the introduction of the new real valued integration variables a^^ and a-^^ , respectively. 
Apart from an irrelevant overall constant, the partition function reduces to 

Z = J P(jPadet^/^((T + ia + m)det^/^((T - la + m) exp ^Tr((T + ia){(x - ia). (5.5) 

As always, the measure of the integral over the matrices a and a is the Haar measure. 
The diagonal mass matrix is denoted by m. This integral greatly simplifies for Nf = 1. 
This case is discussed in Appendix A where we obtain Leutwyler-Smilga sum rules for 
finite values of A^. 



The complex matrix a + ia can be decomposed in 'polar coordinates' as pO 



cr 



+ ia = UAV-\ (5.6) 



where U and V are unitary matrices and A is a diagonal real positive definite matrix. Since 
the r.h.s has Nf more degrees of freedom than the l.h.s., one has to impose constraints 
on the new integration variables. This can be achieved by restricting U to the coset 
U{Nf)/U{l)^f , where [/(l)^-f is the diagonal subgroup of U{Nf). In terms of the new 
variables the partition function reads 

Z = J J{A)VAVUVVdet^^\UAV-^ + m)det^/\VAU~^ + m) exp( ^TrA=^), 

(5.7) 

where the integral over U is over U{Nf)/U{l)^f and the integral over V is over U{Nf). 

For Nf flavors we have Nf condensates which break down the symmetry of the action 
to U{Nf)/U{l)^f leaving us with Nj Goldstone modes. When we allow for fluctuations 
of the topological charge the phase of the determinant becomes massive and the usual 
number of of Nj — 1 Goldstone modes is recovered. 

The term proportional to m plays the role of a small symmetry breaking term. The 
integrals over the nonzero modes will be performed by a saddle point integration at m = 0, 
whereas the integrals over the soft modes will be accounted for exactly for a fixed value 
of mN. The nonzero mode part of the partition function is given by 

Z{m = 0) = y J(A)ciAdet^Aexp( ^TrA^). 

(5.8) 



The leading order contribution in of the integral over A can be obtained by a saddle 
point approximation. The saddle point equations for the A integrals read 



A, = ±A. (5.9) 

The negative solution is not inside the integration manifold and can be omitted. 

At the saddle point in A, the [/—dependence can be absorbed into V. The the 
[/—integration yields a finite irrelevant constant. The remaining integral over V which is 
over U{Nf) can be split into an integral over the phase of the determinant and an integral 
over SU{Nf). We treat m as a small parameter and expand the determinants up to first 
order in m. The result for the m dependent part of the partition function is 

, ' = da PV^expf Tr(m\/"^ exp(-ia/iVf) + m\/exp(m/A^f)), 

Z{m = 0) Jo Jdotv=i 2 

(5.10) 

which coincides with the result for zero topological charge derived by Leutwyler and 
Smilga for the QCD partition function using chiral perturbation theory. As was shown 
in [0 this result can be easily generalized to an arbitrary topological charge sector. 

The value of the chiral condensate can be obtained from eq. (2.6). In the case of equal 
positive masses we have 

< qq >u=o= lim lim < Ti{V + V-') >,=o ■ (5.11) 

For Nf = 1 the integral over V is absent, and the sign of the quark condensate is inde- 
pendent of the sign of m. For more than one fiavor the order of the limits allows us to 
perform the V integral by a saddle point approximation. In the case of equal positive 
masses the saddle point is at \^ = 1 which allows us to identify the parameter S and the 
chiral condensate 

<qq>=E. (5.12) 

This completes the reduction of the partition function. 

As observed in ||^, for more than one fiavor the value of the condensate depends on 
the sign of the quark mass. For example, in the case of two fiavors with equal negative 
masses and ^ = 0, the saddle point is at \^ = —1. 
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6. Spectral density of the random matrix model 



In this section we obtain the spectral density of the random matrix partition function 
(4.1). We use the orthogonal polynomial method developed by Wigner, Dyson and Mehta 
PT| , p2| in the context of the invariant random matrix ensembles. The first step is to 
rewrite the matrix integration in the partition function (4.1) in polar coordinates. For an 



arbitrary complex matrix this amounts to the transformation ^0 



T = UAV 



-1 



(6.1) 



where U and V are unitary matrices and A is a positive definite diagonal matrix. Since the 
r.h.s has n = N/2 more degrees of freedom than the l.h.s., one has to impose constraints 
on the new integration variables. This can be achieved pO| by restricting U to the 
coset U{n)/U{l)"', where f/(l)" is the diagonal subgroup of U{n). The Jacobian of this 
transformation, that depends only on the eigenvalues of A, is given by 



k<l k 



(6.2) 



The integrations over the eigenvalues and the unitary matrices decouple. The latter 
only result in an overall irrelevant constant factor and can be ignored. The eigenvalue 
distribution is thus given by 



Pn(Ai, ■ ■ ■ , A„) = j(A) nn(^^ + ^)) exp(-ns2 y: xd. 

f k 



(6.3) 



k=l 



The spectral density pi(A) is obtained by integration over the remaining n — 1 eigenvalues 

/n 
l[dXkPn{Xi,---,K)- (6.4) 

k=2 

These integrals can be evaluated with the help of the orthogonal polynomial method (see 



2T[| for references). The main ingredient is to write the product over the differences of 



the eigenvalues Vandermonde determinant, i.e. 

1 



k<l 



2\2 



(6.5) 
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which up to a constant can be rewritten in terms of orthogonal polynomials Pk as 



The Pk will be chosen orthogonal according to the weight function 

/ d(A')(A2 + mY' exp(-nE2A')Pfe(A')Pi(A2) = S, 
Jo 

For m — these polynomials are well known, 



(6.6) 



kl- 



(6.7) 



1/2 



(6.8) 



r(iv^ + /, + !); 

where the Lf. ^ are the generalized Laguerre polynomials. 

The determinant can be expanded according its definition. All integrals can be per 
formed immediately by orthogonality and, up to an overall constant, we are left with 



n-l 



pi(A)=2Ev^5: 



(6.9) 



where z is defined by 



(6.10) 



The normalization constant has been chosen such that / dXp{X) = N (remind that N — 
2n). The sum can be evaluated exactly with the Christoffel-Darboux formula (which can 
be found in any treatise on orthogonal polynomials), resulting in 



Pi (A) 



{L:U{z)L:Lt\z) - L^J{z)L:Lt\z)) z^f^'/'eM-z), (6.11) 



r(iV/ + n) 

which, constitutes the exact spectral density of the model (4.1). The microscopic limit is 
obtained by taking N ^ oo while keeping A'^A = x fixed (remember that n — N/2). This 
can be achieved from the asymptotic relation 



hrn i-L^(-)=a;-tJa(2V^), 



(6.12) 



where Ja is the ordinary Bessel function of degree a. One notices that to leading order in n 
the difference in eq. (6.11) cancels. However, using recursion relations for the generalized 
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Laguerre polynomials, this cancellation can be achieved explicitly, and the asymptotic 

relation can be applied to the next to leading order terms. The result for the microscopic 
spectral density is 

Ps{x) = ^-(J^^,(Sx) - JNf+l{^x)JNf-l{^x)). (6.13) 
From the asymptotic relation for the Bessel function 

Jy{z) ~ ( — ) cos(z — —V — —) for z ^ oo (6.14) 

we find that 



5] 

lim psix) = — (6.15) 

X— »oo TV 

which agrees with the Banks-Casher relation. As follows from the small x behaviour of 

^^^"^^ iV,!(iV, + l)! ^'-^'^ 

small eigenvalues are strongly suppressed for an increasing number of flavors. 

The formula (6.13) reproduces all diagonal sum rules of Leutwyler and Smilga, e.g., 
the sum 

can be converted into an integral over the microscopic variable x = XN resulting in 

psix)dx _ Ti2p-l)TiNj-p+l) 

x^P \2j T{p)T{p+l)T{Nf+p) ^ ' 

The above spectral density thus summarizes all sum rules and yields new sum rules, e.g. 
sum rules for noninteger values of p. 

It is possible to derive all higher order spectral correlation functions which generate an 
infinite family of new spectral sum rules. An explicit expression for the level correlation 
function is given in [Q , and it generates sum rules obtained previously by Leutwyler and 
Smilga [0]. 
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7. Spectral density for the instanton liquid 



The sum over field configurations in the partition function (2.2) can be approximated 
semiclassically by that of a liquid of instantons. Instead of averaging over all gauge field 
configurations, we average over the collective coordinates of the instantons only, whereas 
1-loop quantum fiuctuations about the instantons are included in the measure. The action 
in (2.2) is the instanton action, which also includes the interaction between instantons. 



We use the so called streamline p4|, |25|| interaction supplemented by a core in order to 
stabilize the instanton liquid. The fermion determinant is calculated in the space spanned 
by the fermionic zero modes with overlap matrix elements that can be derived from the 



streamline configuration |15|. More details on the above instanton liquid model can be 



found in [27 



The numerical simulations were carried out for a liquid of 64 instantons in a Euclidean 
space time volume of (2.378)^ x 4.756 in units of Ag^^. Averages were obtained from 
1,000 statistically independent configurations for Nf = 1, 2, 3 and 100,000 for Nf = 0. 
Our main results are presented in Figs. 1 and 2. The number of fiavors is shown in the 
label of the figures. In Fig. 1 we show the average spectral density n(A) of the Dirac 
operator in the space of the zero modes (the normalization is n{X)dX = 1, and A is 
in units of Aqcd, so ^(A) differs from < p(A) > by a normalization factor). Note that 
for Nf = 2 the thermodynamic limit of the slope n'(0) is zero, whereas for Nf = 3 it is 



negative. This agrees with a recent result of Smilga and Stern [26| according to which the 
flavor dependence of n'(0) is given by ~ {Nj — A)/Nf. 

For Nf = 0, 1, 2 we observe a clear separation of the microscopic scale near zero 
virtuality and the overall average spectral density. The latter is obtained by a flt with a 
smooth function ignoring the region near zero virtuality. This allows us of obtain a flnite 
value for < p(0) >. The microscopic limit of the spectral density is shown in Fig. 2. The 
full line is the spectral density of the the unfolded spectrum < p{fi) > (see ref for an 
exact deflnition), for which the average spacing between the eigenvalues is equal to one. 
The dashed curve shows the asymptotic result for the microscopic spectral density in the 
same units which agrees well with the random matrix result. This is even more remarkable 
if one realizes that the average spectral density is very different from the random matrix 
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result (a semi-circle). It should be stressed that the microscopic spectral density contains 
no free parameters. For Nf = 3 it is no longer possible to extract the thermodynamic 
limit of < p(0) > and an unambiguous comparison with the microscopic spectral density 
cannot be made. The strong supression of small eigenvalues is consistent with eq. (6.16). 

We observe that the spectrum is very stiff. The peak occurs because the smallest 
eigenvalue is always close to its average position. As is also the case for the other invariant 
random matrix ensembles, the fluctuations are much less than for a random sequence of 
eigenvalues, although it looks like that the oscillations in the asymptotic spectral density 
are not reproduced by the liquid of instantons. This may be a finite size effect. 

The Leutwyler-Smilga sum rules can be checked for a liquid of instantons. The results 
completely agree with the the exact results in 0] (see for a detailed comparison). 

8. Conclusions 

We have studied the spectrum of the QCD Dirac operator near zero virtuality with the 
help of the microscopic spectral density, which is obtained by enlarging the spectrum by 
V4 as V4 ^ oo. The existence of this limit is an immediate consequence of the spontaneous 
breaking of chiral symmetry. 

We have conjectured that the microscopic spectral density is universal. This should be 
contrasted with the average spectral density which depends on the details of the system. 
The latter varies on a scale of 0{V^), whereas the first shows fluctuations on the micro- 
scopic scale of C(l/V4). The implicit assumption is that both scales can be separated. 
On the basis of the universality we have formulated a random matrix theory that apart 
form the chiral structure of the QCD Dirac operator has no other information content. 

It was shown that this random matrix theory is equivalent to the extreme low-energy 
limit of the QCD partition function. Therefore all Leutwyler-Smilga sum rules are re- 
produced. Using the orthogonal polynomial technique developed by Wigner, Dyson and 
Mehta, we were able to obtain the exact analytical result for the microscopic spectral 
density. 

In general it is very hard to obtain exact numerical spectra for a field theory. However, 
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in the instanton liquid model for the QCD vacuum we were able to obtain sufficient 
statistics for the microscopic spectral density, and we found a very satisfactory agreement 
with the random matrix model. In view of this, it would be very interesting to compare 
our results to spectra obtained from lattice QCD. 



Appendix A 

In this appendix we derive the sum rules for one flavor and finite values of N. For 
Nf — 1 the partition function (5.5) simplifies to 

, nj:^ 

dada[a + ia — m) ^ [a — ia — m) exp( — (cr + ia){a — icr)), (A.l) 

where an irrelevant overall constant has been suppressed. At finite N the pre-exponential 
factors can be expanded as a binomial series which provides us with an expansion in 
powers of m. The coefficients are elementary integrals, and for the m— dependent part of 
the partition function we find 

It should be noted that no approximations have been made. On the other hand, the 
fermion determinant can be written as a product over the eigenvalues which leads to the 
expansion 

where the average is with respect to the massless partition function. By equating the co- 
efficients of the powers of we obtain sum rules for the inverse powers of the eigenvalues 
that are valid for any value of N. The sum rule (3.3) for Nf = 1 is reproduced. It is valid 
for any value of A^. The second sum rule is modified by the factor (1 — 2/N). For N — 2 
we find zero which is correct because in this case there are no terms that contribute to 
the sum n' ^ n in (A. 3). 
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Figure Captions 

Fig. 1. The eigenvalue density n{X) for Nf — 0, 1, 2, and 3. The area below each curve 
is normalized to 1, and the binsize is 0.00125 for Nf = and 0.0025 in the other cases. 

Fig. 2. The microscopic hmit of the spectral density in Fig. 1. We show the density 
for the unfolded spectrum where the average position of the n'th eigenvalue is at = n. 
The dashed line represents the exact microscopic spectral density for the same 

normalization, i.e. E = tt. 
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